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The  particle  trapping  phenomenon  that  occurs  in  the  Whistler 
emission  process  is  more  complex  than  the  trapping  by  an 
infinite  wave  in  a  homogeneous  system.  Because  of  the 
inhomogeneous  background  magnetic  field  and  the  antiparallel 
electron  and  wave  velocities,  electrons  experience  trapping  for 
only  a  finite  length  of  their  orbits.  Successive  trapping  and  - 
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de trapping  snrsots  occur ,  avan  for  an  infinite  wave  train.  A 
aelf -cone latent  theory  for  the  emission  process  requires 
knowledge  of  the  orbits  throughout  these  different  phases.  One 
snist  follow  the  electrons  for  many  trapping  periods  and 
accurately  track  the  phase  so  that  the  detrapping  point  can  be 
computed .  The  present  report  describes  an  asymptotic  theory 
that  allows  this  to  be  done,  basically  by  ^finding  an  adiabatic 
invariant  of  the  trapped  electron  motion. ^The  calculation  is 
done  by  an  asymptotic  ordering  which  relatjes  the  various  small 
parameters  in  order  to  bring  the  appropriate  physical  processes 
into  the  expansion  in  a  workable  way.  It /then  turns  out  that 
this  ordering  corresponds  very  well  to  the  relations  that  occur 
in  practice,  with  the  numerical  value  of  /the  basic  expansion 
parameter  being  about  .1  .  This  suggests  that  the  analytic 
theory  will  have  quite  good  accuracy  and /be  useful  for  practical 
applications.  / 
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ABSTRACT 


The  particle  trapping  phenomenon  that  occurs  in  the  Whistler  emission 
process  is  more  complex  than  the  trapping  by  an  infinite  wave  in  a  homogen¬ 
eous  system.  Because  of  the  inhomogeneous  background  magnetic  field  and  the 
antiparallel  electron  and  wave  velocities,  electrons  experience  trapping  for 
only  a  finite  length  of  their  orbits.  Successive  trapping  and  detrapping 
events  occur,  even  for  an  infinite  wave  train.  A  self-consistent  theory  for 
the  emission  process  requires  knowledge  of  the  orbits  throughout  these  dif¬ 
ferent  phases.  One  must  follow  the  electrons  for  many  trapping  periods  and 
accurately  track  the  phase  so  that  the  detrapping  point  can  be  computed. 

The  present  report  describes  an  asymptotic  theory  that  allows  this  to  be 
done,  basically  by  finding  an  adiabatic  invariant  of  the  trapped  electron 
motion.  The  calculation  is  done  by  an  asymptotic  ordering  which  relates 
the  various  small  parameters  in  order  to  bring  the  appropriate  physical  pro¬ 
cesses  into  the  expansion  in  a  workable  way.  It  then  turns  out  that  this 
ordering  corresponds  very  well  to  the  relations  that  occur  in  practice,  with 
the  numerical  value  of  the  basic  expansion  parameter  being  about  .1  .  This 
suggests  that  the  analytic  theory  will  have  quite  good  accuracy  and  be  use¬ 
ful  for  practical  applications. 


I.  INTRODUCTION 


The  report  describes  so  asymptotic  orbit  calculation  for  high  energy 
electrons  moving  along  the  earth's  magnetic  field  in  the  presence  of  a 
large  amplitude  Whistler  wave.  This  calculation  is  the  necessary  first 
step  in  the  development  of  a  self-consistent  non-linear  theory  of  the 
triggered  emission  process.  An  adiabatic  invariant  for  the  trapped 
particle  motion  is  found  from  a  multiple  time  scale  analysis  of  the 
trajectories.  This  allows  the  orbits  to  be  followed  for  many  trapping 
periods  to  keep  track  of  the  phase  and  to  compute  the  detrapping  point. 
From  this  an  expression  for  the  nonlinear  resonant  current  can  be  written 
down  to,  in  principal,  close  the  WKB  wave  equations  describing  the  emis¬ 
sion  and  propagation  characteristics. 

In  section  II,  the  orbit  theory  is  developed  as  an  asymptotic  expan¬ 
sion  to  bring  out  the  main  physical  characteristics.  The  resulting  order 
ing  is  compared  with  typical  experimental  parameters  in  section  III.  Sec 
tion  IV  reviews  the  properties  of  the  orbital  phase  space,  including  an 
evaluation  of  the  trapped  particle  separatrlx. 


II.  ORBIT  ANALYSIS 


In  the  cantered  dipole  approximation  (sea  Helllwell's  book  Whistler 
and  Related  Jonospharic  Phenomena,  Appendix)  the  magnetic  field  strength  is 
given  by 
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where  $  is  the  geomagnetic  latitude  and  R  the  geocentric  radius  which  can 
be  expressed  as 
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In  the  region  of  interest  about  the  equator,  i.e.  for  small  $,the  field  can 
be  expressed  as 


B  -  B  (1  +  *- 
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s  being  the  arclength  along  the  line  and  L  being  the  scale  length  of  the 
gradient.  In  this  field  the  electron  equations  of  motion  in  the  presence 
of  the  Whistler  field  are,  [1] 
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B  stands  for  the  wave  magnetic  field  amplitude. 


Where  the  electron  gyrofrequency  is  given  by 


(8) 


the  magnetic  moment  y  is  defined  as 


(9) 


where  is  the  magnitude  of  the  electronvelocity  vector  perpendicular  to 
the  earthmagnetic  field  direction.  The  phase  $  of  this  vector  has  been 
replaced  by  the  relative  phase  y  between  v^  and  the  electric  field  vector 
of  the  triggering  wave: 
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(10) 


In  order  to  arrive  at  this  set  of  differential  equations  for  the  behavior 
of  an  electron  we  have  utilized  the  smallness  of  the  gyroradlus  compared 
to  the  scale  length  of  the  magnetic  field. 

This  makes  y  in  the  absence  of  the  wave  a  constant  of  the  motion. 

In  the  presence  of  the  wave  both  u  and  the  total  particle  energy  will 
change.  Adding  equations  (5)  and  (6)  gives  the  exact  expression  for  the 
change  of  total  kinetic  energy: 
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Rigorously  speaking  everything  on  the  right-hand  side  of  equation  (II) , 


except  for  ft  and  B  ,  is  a  function  of  s  and  t  .  The  evaluation  of 
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vx  is  described  by  equation  (6)  which  can  be  written  alternatively,  as 


(12) 


Furthermore,  the  phase  velocity,  ^  ,  of  the  wave  will  change  as  a  function 
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of  s  and  t  and  after  tha  instability  touches  off  B  ,  the  wave  field,  will 
grow  in  space  sod  tiaa. 


While  these  changes  occur  on  a  slow  tine  scale  the  angle  Y  will 
change  very  rapidly  for  aost  electrons.  Since  the  propagation  character¬ 
istics  of  the  waves  are  linear  (due  to  cold  electrons)  and  the  wavelength 
and  period  of  oscillation  are  saall  cospared  to  the  spatial  and  temporal 
scales  of  evolution,  the  Whistler  wave  can  be  described  by  a  WKB  or 
eikonal  form, 

-  !♦„(*. t) 

E  -  P  E(s,t)e  *  (13) 

where  ?  is  the  polarization  vector,  E  the  wave  amplitude  and  $w  the  wave 
phase.  The  last  is  determined  from  an  Integration  of  the  ray  equations 


(14) 


Then  the  phase  Y  between  the  wave  and  the  particles  becomes 

T  *  ♦„(»,«)  +  (15) 

where  the  particles  phase  $  is  evolving  according  to  the  equations  of  motion: 


the  rate  of  change  of  Y  is  given  by 

dY  ds  3*w  ,  3<*w  (  d» 
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The  homogeneous  field  calculation  with  its  resonance  velocity 

V  -  iqJ  Helliwell's  experiments  show  that  relevant  wave  phase  velocities 
rM  k 

~  are  of  the  order  of  the  parallel  velocity  of  the  group  of  electrons  of  inter¬ 
est.  Of  course  the  plasma  under  consideration  is  trapped  in  the  mirror  magnetic 
geometry  of  the  earth's  dipole  field  which  lends  itself  to  loss  cone  driven  in¬ 
stabilities.  However,  observations  indicate  that  the  medium  Is  only  weakly 
unstable,  linearly.  This  leads  us  to  conclude  that  even  though  an  inbalance 
between  and  vx  exists  both  can  still  be  on  the  same  order  of  magnitude. 
Putting  all  of  this  together  says  chat  the  first  term  on  the  right  of  equation 
(7)  will  be  of  order  e2  compared  to  the  rest.  And  we  can  write 

— ^  ■*  kv  —  bi  +fl  +  0(e2)  (21) 

dt 


This  says  that  for  certain  v^  the  particles  will  not  change  their  phase  rela¬ 
tion  with  respect  to  the  wave  field.  This  occurs  for  velocities  near  the 


resonant  velocity 


V  -  ~  (22) 

res  k 


It  is  convenient  to  define  a  new  variable 


v.  «  V  +  v  with  —  <<  1  (23) 
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and  one  can  write 
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The  phenome non  of  Cropping  is  contained  in  equation  (5)  when  the 
lnhomogenait#  4 
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B  —  can  becoae  sufficiently  small  near 
0  L2 

BW 

the  equator,  s  -  0  ,  that  the  wave  term  -  ^  —  vA  cos  Y  dominates  and 

o 


causes  the  electrons  to  oscillate  about  some  stable  phase  at  a  frequency 

,w  v  I/2 


"tr-  (°o  ^7  ) 
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the  so-called  trapping  frequency  u>tf  for  well  trapped  particles. 

We  now  focus  attention  on  the  trapped  particles,  and  develop  an 
asymptotic  ordering  of  parameters  that  permits  an  expansion  of  the  electron 
orbits.  We  will  show  later  that  this  ordering  accurately  describes  the 
real  parameters  in  the  observed  emission  process.  Re-writing  in  terms 

of  the  equatorial  gyroradius  p  ■  — 
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tr 
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B 

one  notes  first  that  u  is  much  smaller  than  0  since  —  is  a  very  small 

tr  o  d 

_  o 

number  (in  practice  on  the  order  of  10  5).  In  practice  kp  is  order  1  and 
thus  should  remain  unordered.  Then  we  can  define  a  smallness  parameter  c 
as 


tr 
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indicating  at  the 


V 

time  that  g-  ~  e2 
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Since  Y  is  intrinsically  an  order  one  quantity  we  find  for  the  region 
of  Interest  (jr  sufficiently  small  so  that  both  terms  on  the  right-hand  side 

.  ^ij( 

of  equation  (5)  can  compete)  gg  -  u>tr  Y  .  Note  that  for  untrapped,  nonres¬ 
onant  particles  the  phase,  Y,  changes  on  the  gyroperiod  time  scale  and  another 
expansion  procedure  would  have  to  be  used. 


In  terms  of  this  now  variable  the  equations  of  notion  read 
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where  from  the  definition  of  V  — « 
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was  used. 
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For  trapping  to  occur,  we  must  have  the  Inhomogeneity  and  wave  terms  in  equa¬ 
tion  (26)  comparable.  The  terms  on  the  right-hand  side  of  this  equation  scale 
successively  like  this: 


a  B»  v.2s  fl»  * 
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When 


treating  a  constant  frequency  wave  train  the  last  term  ^  can  be  neglect¬ 


ed  since  V  is  not  an  explicit  function  of  time.  The  proportionality  (30) 


can  be  re-expressed  as: 
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By  assumption  —  «1  such  that  the  last  term  in  (31)  is  very  small  compared 
v* 

to  its  precessor  — -  *>  1  was  used  again.  Also  using  kp  *  1  again  it  is  true 
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that - * - .  Then  (31)  demands  essentially 
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(32) 


Now,  we  also  wane  the  a  notion  to  bo  on  a  slower  tine  scale  than  the 
bouncing  in  the  wove  crests  so  that  s  is  essentially  forzen  (adiabatically) 
while  trapping  notion  ensues.  Therefore 


—  -  e°  utr  where  0  <  a  <_  1 
s 

:  v,  &  i  l 

has  to  be  determined  self  consistently.  —  scales  like  —  ~  - 
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or 
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From  both  (32)  and  (34)  one  can  solve  for 
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(35) 


In  order  to  be  able  to  treat  B  -  Bq  on  both  these  scales  as  essentially  con- 
stant  we  want  the  first  correction  tern  ^  — ^  to  cone  in  on  the  next  slower 
time  scale  (ea)  .  Hence 


(36) 


Both  (35)  and  (36)  can  be  solved  to  give 
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(37) 


This  completes  the  asymptotic  ordering,  using  the  small  parameter  c  ■  (Bw/B  ) 

o 

The  fast  tine  scale  for  trapped  particle  is  u  -1  on  which  scale  s  is 

1/3  ~l 

frozen.  On  this  and  the  next  slower  time  scale  e  u>  ,  where  s  motion 
occurs,  the  magnetic  field  can  be  treated  as  constant.  It  follows  from  (32) 


and  sine*  v  changes  as  t  on  the  fast  tins  scale  w 


n  v  a  v.  a  v, 

o  i-  O  1 


v 

—  -  e 


v 


i 


(39) 


Thus  the  corrections  to  v.  -  V _  are  coming  in  on  third  order  in  the 

»  res 

l  /3 

expansion  parameter  e  ,  and  are  not  needed  for  computing  the  motion  along 
the  magnetic  lines.  To  lowest  order  the  equations  of  motion  for  trapped 
particles  become: 
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In  proceeding  toward  a  solution  of  this  set  of  equations  it  can  be 
noted  first  that  the  time  derivative  of  equation  (43)  becomes 
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The  last  term  in  this  equation  scales  like 
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Comparing  this  to  the  left-hand  side  which  scales  like 


one  finds 
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thee  this  term  Is  s  third  order  correction  in  c  leaving 
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One  recognizes  the  first  term  on  the  right-hand  side  as  the  wave  ten, 

w2  cosY  ,  and  the  others  as  inhomogeneity  terms.  We  define  an  Inhomogeneity 
c* 

factor,  S  ,  according  to. 
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Hence,  finally  the  Y -equation  takes  the  form 

— -  +  w2  cosY  ■  u>2  S 

tr  tr 


There  exist  bounded  oscillatory  solutions  to  this  equation  as  long  as  the 
inhomogeneity  on  the  right-hand  side  satisfies 

{S}  <  1  (52) 


assuming  for  now  that  the  frequency  stays  constant. 


an  oscillator  aquation  with  oscillation  frequency 
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Thus  Y,  becomes 
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Since  4  Is  Inherently  secular  on  the  fast  tine  scale  one  has  to  compute 

^  to  requisite  order  prior  to  being  able  to  evaluate  4  •  To  find  ■— &- 

the  next  equation  in  the  expansion  hierarchy  is  needed.  It  reads 
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Our  objective  is  to  eliminate  secularlties  of  Y.  on  the  fast  time- 
-i  32po 

scale  «>o  .  Based  on  equation  (58)  ■■  vanishes  eliminating  the  second 

term  on  the  right-hand  side  of  equation  (65).  Also  the  third  and  fifth  term 
cannot  drive  the  left-hand  oscillator  at  frequency  wq  .  The  fourth  term: 
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is  driving  at  twice  the  frequency  u>0,  also  dropping  out  of  interest.  It  is 
then  left  to  require  that 
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will  not  drlvs  st  <i»o  lending  to  s  fsst  tins  scsls  ssculsrity  in  T2  .  This 
will  b«  gusrsntssd  under  ths  following  two  conditions: 


snd 


Au  •  0 
o 


(68) 


(69) 


Therefore  the  combination  Awq  ,  the  amplitude  of  trapping  oscillations  times 
the  frequency  is  constant  on  the  slow  time  scale.  This  is  a  rather  unexpected 
constant  of  the  motion  pertaining  to  the  slow  time  evaluation  of  the  trapped 
particle  motion.  And  an  additional  fact  is  of  great  importance:  equation  (69), 
which  allows  us  to  Integrate  A  as  follows 


A 


dt  «tr  /  sin[-po(s(t))] 


(70) 


From  the  equations  (AO),  (SO),  snd  (60)  everything  is  known  to  compute  A  • 
Assembling  these  results  we  found  an  analytic  description  of  well  trapped 
electrons. 


f  -  PQ  +  A  sin  A 


A 


dt  w  /sin(-p) 
tr  o 


(71) 

(72) 


With  the  adiabstlc  invariant 


A  u  ■ 


constant 


(73) 


From  aquation  (43)  vo  find  that  tha  par  cur  bad  velocity  v  Is  of  constant 
amplitude 

**\> 

V  ■  -  cost 

k 

Summarizing  thaaa  results  tha  trappad  particles  are  found  to  gyrate 
In  vortices  in  v  -  f  -  space  on  the  fast  tlsescale  .  Excursions  in  v 
are  of  constant  amplitude.  The  average  phase  angle  PQ  is  slowly  changing 
as  tha  particles  are  drifting  further  along  the  field  lines.  Outside  the 
region  of  trapping  |S|__>  1  »  the  particles  follow  adiabatic  orbits  in  an  in¬ 
homogeneous  field  to  leading  order.  For  simplicity  an  infinite  wave  train 
of  constant  amplitude  was  considered. 
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III.  COMPARISON  WITH  EXPERIMENTAL  PARAMETERS 


After  bmvimg  developed  the  analytic  expression*  for  trapped  electron 
orbits  the  question  erlses,  how  do  our  assumptions  compare  with  the  measure¬ 
ments?  Are  the  expansion  paraaeters  apt  to  validate  the  calculations? 

In  order  to  check  the  expansion  against  the  data  given  in  the  litera¬ 
ture  our  interest  will  focus  on  a  field  line  designated  as  L  «  4  . 

The  L  number  is  the  multiple  of  earth  radii  (Rq  *  6371  km)  that  gives 
the  straight  distance  from  the  center  of  the  earth  to  the  equator  of  the 
fieldline. 


Granted  a  dipole  representation  of  the  earth  magnetic  field  the  1*4 
line  originates  at  ■  60°  latitude  on  the  earth's  surface.  This  corresponds 
to  a  fieldline  length  of 


-  57,360  km  (74) 

(The  arc  length  in  this  form  was  given  by  Chapman  and  Sugiura  in  1956)  [2]. 
At  and  in  the  vicinity  of  the  equator  of  the  fieldline  under  consideration, 
the  magnetic  field  can  be  expressed  as 

,',o(1+^)  •  <75) 


The  scale  length  is  given  by 

/2”L  R 

L  -  - -  -  12,013  km  (76) 

3 

For  the  proposed  comparison  of  scales  one  needs  to  compute  the  relevant  plasma 
parameters  like  fl  ,  u)  ,  V  ,  p  ,  etc.  as  well  as  wave  parameters. 


The  dipole  approximation  to  tho  earth  aagnatlc  flald  givas  an  alactron  gyro- 
frequency  of  ■  13.63  kHs  at  an  altltoda  of  3  aarth  radii  above  tha  earth's 
•qua tor  [2].  Thia  masher  compares  wall  with  tha  more  datailad  analysis  by 
Fougara  using  spharieal  harmonics  which  givas  f  ■  14  kHz  [3].  Tha  corres- 
ponding  “  2ir  f  ■  88000  rad/ssc.  - 

Typical  obsarvation  are  parfornad  with  tha  fioldstrangth  of  tha  trigger¬ 
ing  signal  lying  in  tha  range  of  1  -  10  ay  [4].  The  commonly  used  unit  of  y 
is  defined  by  ly  -  10* 9  Tesla  ■  10~5  Gauss. 

Hence: 

£ .  s  ■up.IIx  .  l0-s . 

Bq  500y 

This  proves  to  be  a  very  good  expansion  parameter.  It  is  known  experimentally 
that  the  frequency,  where  triggering  of  emission  takes  place  most  frequently ,  is 

t 

half  the  minimum  gyrofrequency  [5]  of  electrons  along  the  fleldllne  0q  : 

n 

u >  z.  —  *  44000  rad/sec;  f  *  7  kHz 
2 


Helliwell's  measurements  of  the  ',nosa"-frequcncies  of  risers  and  fallers  allow 
to  conclude  on  the  plasma  frequency  [2]  which  he  gives  as 


u>pe  »  7.6  x  10  rad/sec  corresponding  to 


180  electron/cm3  .  How  we  are  in  the  position  to  compute  the  wavelength  based 
on  the  linear  dispersion  relation  for  whistler  mode  signals  in  homogeneous 


magnetized  plasmas: 


c  -  1  - 


2  2 
c  k 

u.2 


u>  2 


u r 


-  0 


(77) 


u 


1  - 


ce 


Hi 


where  w  is  the  electron  gyrofrequency.  Plugging  in  the  numbers  k  turns  out 


k  »  2.Jr-  corresponding 


2w 


to  s  wave  length  X  of  approximately  X  »  —  »  2.5  ka  .  We  also  take  the 

k 


resonance  velocity  on  the  order  of  the  particle  velocities 


w-fl 


^re.—"  ^  •  ’  Vvx 


(78) 


(This  is  also  justified  by  a  set  of  nodal  parameters  in  a  late  paper  of  R.A. 
Hellivell  and  U.S.  Inan  [6].  Fron  this  the  gyroradius  p  is  0.2  km.  a  tenth 
of  the  wave  lenth.  Then,  very  closely,  kp.  appears  to  be  an  order  1  quantity. 


The  definition  of  the  crapping  frequency  was 


-  n2  Ko  —  -  200  rfr 


rad 


tr 


sec 


To  within  the  accuracy  that  kp  is  1  it  will  be  found  that  “tr/®0  *  e  •  Here 


e  is  being  defined  by  the  square  root  of  the  ratio  of  the  field  strengths. 


si  B 

Next  wc  are  to  compare  the  quantity  ^  ££  with  —  to  find  an  upper 

o 


bound.  Ue  can  take  the  length  Lp  of  the  particle  interaction  region  FIR  as 


estimated  by  Helliwell  [7] 


L  -  800  km 
P  — 


or  the  arc  length  s  up  to  which  the  effect  of  particle  trapping  can  be  found. 
I.E.  s  for  which  the  two  terms  on  the  right-hand  side  of  the  equation  of  motion 


for  v  balance.  The  condition  is 

■ 


It 

f  2V« 

- k  ’’Ml. 

J 

^  Lw2 
tr 

L  '' 

(79) 


from  which  follows  that  s  can  range  up  to  1500  km  . 


With  Helliwell' s  estimate  we  find 


-  2.2  x  10 


nW  _5 

which  it  amasingly  dose  to  —  *  10  ,  verifying  our  assumption 

B 

o 

si  B 

- -  —  '  e2 

LkL  Bq 

Further  we  set  —  —  -  e°  <o  .  By  definition  —  ■  v,  .  There  follows 
s  dt  cr  dt 


21.7  -  • 

a  dt  s  9 

For  our  purposes  of  finding  analytic  expressions  for  the  orbits,  we  expand 

1/3  l 

hence  in  terms  of  e  "xo  *  To  conclude  this  section  we  give  a  short 
table  of  parameters  relevant  to  the  orbit  calculations: 


weq  -  760,000  rad/sec 
pe 


L  *  12,013  km 


•  88,000  rad/sec 


X  -  2.5  km 


44,000  rad/sec 


k  -  2.5  1/km 


200  rad/sec 


0.2  km 


-  2.2  x  10 


V»  '  Vx  '  vres“  17,440  ka/* 


e,/3  -  1/10 


1  -  4 


1  -  10  m  y 


R  -  6,371  km 
o 


0  *  1500  km 


IV.  SEPARATE!!  ANALYSIS 


With  the  solutions  to  ths  equations  of  notion  the  wave-particle  problem 
has  been  solved  on  a  microscopic  level.  We  found  that  the  Individual  electron 

A 

executes  an  orbit  which  is  parametrized  by  two  quantitites  A  and  $.  A  is  the 

A 

amplitude  of  trapping  oscillation  in  Y  and  A  the  initial  phase  on  this  orbit, 
as  it  appears  in  the  Integration  of  #  (eq.  70  ). 

Giving  a  pair  of  these  parameters  allows  one  to  integrate  the  orbit  back 
or  forward  to  the  entrapping  or  detrapping  point.  These  points  are  determined 
by  the  intersection  of  the  separatrix  in  v  -  Y  -  space  with  the  individual 
orbit.  Exterior  to  those  two  points  the  electrons  follow  adiabatic  orbits  in 
the  geomagnetic  field  to  lowest  order. 

At  exactly  the  equator,  S  ■  0  ,  the  medium  appears  homogeneous  and  the 
phase  space  picture  is  as  shown  in  Figure  1. 
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Figure  1.  , 

Trajectories  of  cyclotron  resonant  electrons  in  the  v-T  plane  at  the  equator, 

S-0,  or  for  the  case  of  a  homogeneous  medium.  The  open-ended  lines  correspond 
to  untrapped  particles  while  the  closed  trajectories  represent  trapped  particles. 
The  separatrix  is  labeled  by  a-1  . 


Clearly  the  closed  line  orbits  of  trapped  electrons  are  separated  from 

r~  “tr 

the  open-ended  lines  for  v  greater  than  /2  — —  . 

k 

The  equation  for  these  orbit  paths  is 


1  k2v2 


+  sin  Y  -  a 


tr 


(82) 


Each  path  is  labeled  by  a  different  constant  of  integration,  a,  of  the  first 
(energy)  Integral  on  the  right-hand  side  of  equation  (82) .  For  this  constant 


being  equal  to  1  corresponds  to  the  equation  for  the  separatrix. 
Away  from  the  equator,  the  equation  for  the  orbits  is 


Bounded  orbits  sppssr  for  0  <_  |s|  £  1  .  Again  the  choice  for  the  constent 
of  integration,  e,  specifies  each  trajectory.  A  number  of  them  ere  depicted 
in  Figure  2.  for  S  •  0.5  . 


Figure  2. 

Trajectories  of  cyclotron  resonant  electrons  in  the  v-f  plane  for  S-0. 5  . 

The  separatrix  of  the  closed  trajectories  about  PQ  for  -w/ 2  <  PQ  <  0  would 
be  parametrised  by  aa#p-0.342  . 

The  constant  of  integration,  a.  and  its  functional  dependence  on  S  can 
be  found  graphically  and  analytically  after  setting  v  *  0  .  This  corresponds 
to  the  intersections  of  the  contour  lines  with  the  Y-axis.  The  intersection 
values  of  V,  are  the  solutions  to  the  transcendental  equation 

a  +  St  -  sint  (84) 

The  straight  line  with  slope  S  on  the  left  equals  the  sine  function  on  the 
right  for  several  y's  depending  on  the  choice  of  a  .  This  is  illustrated  by 
Figure  3. 
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Figure  3. 

Graphical  determination  of  the  constant  of  integration  for  the  separatrix 
and  its  width  in  ¥  .  See  text  for  details. 

Since  ¥  is  the  angle  between  vx  and  the  electric  field  vector  of  the 
wave  it  is  sufficient  to  look  at  one  interval  of  width  2tt  on  the  ¥  axis. 
For  0  <  S  <_  1  there  will  be  maximally  3  intersections  for  small  a's,  two 
of  which  define  the  width  in  ¥  of  the  corresponding  trapped  particle  orbit 
the  third  one  belongs  to  an  open  trajectory.  The  condition  for  the  separa 
trix  is  that  the  open  trajectory  and  the  closed  trajectory  just  touch.  Th 
happens  for  a2  where  the  straight  line  a^  +  S¥  becomes  tangential  to  sin  ¥ 
(see  Figure  3).  The  value  ¥^  where  this  happens  is 

¥_  -  cos"1 S 


with  0  <f  <.  w/2  for  S  >  0  .  The  analog  holds  for  S  <  0  .  Tha  value  of 
a  correspond lag  to  tha  saparatrlx  is  givan  by 

aMo  -  S  coa-1S  -  A  -  S2  < 

sap 

Hance,  the  equation  for  the  separatrlx  about  the  trapped  particle 
island  cantered  about  ,  with  -ti/2  _<  PQ  <,  0  ,  is 


S  cos”1  S  -A  -  S2  +  S4»  -  sinf 


This  allows  us  to  compute  the  entrapping  point  for  nearly  resonant  electrons 
streaming  into  an  infinite  constant  amplitude  triggering  wave.  The  equation 
to  be  evaluated  is  obtained  by  inserting 


v  «  (AuQ)  co 


u»tr  A  -  Sz(t) 


*(J  “tr  n  -  S(t)  dt)  . 

— -  ,ln([  “tr  *  S(t)  4t) 


into  equation  (85)  and  finding  the  root. 

The.  separatrlx  has  roughly  the  shape  of  an  ellipse  as  can  be  seen  from 
Figure  2.  In  order  to  find  an  analytic  expression  for  its  width  in  T  as  a 
function  of  S  it  is  necessary  to  reduce  the  transcendental  character  of  the 
sine-function  involved  in  equation  (84). 

It  was  found  that  a  good  approximation  for  sin  4*L  is  given  by 

/  12 

-  —  Y? - y  -  8  ■  a  -  S4L  (90) 

?  L  L  sep  »» 


in  the  interval - <4*.  <  -  rr  and 

2  ~  L 


25  - 


(91) 


_4 

IT 


£  +  ii.  m  m  m  -  ST 
Z  •  L  L  sep  S  L 


In  the  interval  -*  _<  T^  <_  0  .  The  solution  to  either  one  of  these  equations 

gives  the  left-hand  intersect  of  the  seperetrlx  contour  with  the  T  axis.  The 
right-hand  intersect  is  given  by  equation  (85).  Therefore  the  width  of  the 
seperetrlx  ATggp  ■  TR  -  TL  is  given  by 


AT  -  cos-1  S  + 
sep 


2  8  a  '  „2 


(92) 


A  +  7177“) 

2  8  8  *  it2  \2  *•**/ 


The  upper  expression  Is  valid  for  the  Interval  -  —  <  T.  <  -  ir  ,  the  lower  ' 

2  —  l  — 

expression  for  -ir  <_  <_  0  .  The  width  in  v  of  the  separatrix  Av  is  given  by 

Av  -  2  —  /S  cos"1  S  -  /l  -  Sz  (93) 

k 

So  far  S  was  treated  as  a  positive  quantity.  It  is  straight  forward  to  write 
down  the  expressions  for  -1  <  S  <  0  . 

Belabouring  the  separatrix  to  this  extent  is  of  importance.  All  of  the 
interior  are  trapped  particles  which  represent  the  dominant  contribution  to 
the  resonant  current  density.  This  current  in  turn  is  responsible  for  the 
modifications  of  the  triggering  wave,  the  emission.  The  next  step  to  be 
taken  is  to  Integrate  over  the  resonant  phase  space  region  to  compute  this 
current.  It  will  be  to  our  advantage  that  the  wave  amplitude  is  growing  on 
a  slower  time  scale  than  the  trapping  motion. 
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